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We consider two-dimensional Ising models with randomly distributed ferromagnetic bonds and 
study the local critical behavior at defect lines by extensive Monte Carlo simulations. Both for 
ladder and chain type defects, non-universal critical behavior is observed: the critical exponent of 
the defect magnetization is found to be a continuous function of the strength of the defect coupling. 
Analyzing corresponding stability conditions, we obtain new evidence that the critical exponent f 
of the bulk correlation length of the random Ising model does not depend on dilution, i.e. = 1. 
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The presence of quenched randomness may drastically 
change the critical properties of magnetic systems. For 
disorder which is coupled to the energy density, i.e. in 
particular for random bond and random site dilution, 
the relevance-irrelevance of the perturbation at a second- 
order phase trajisition point is given by the well known 
Harris criteriond. If the specific heat exponent of the 
pure system is positive, a > 0, a new random fixed point 
is expected to control the critical properties of the di- 
lute model. The marginal situation in the Harris crite- 
rion, a = 0, is represented by the two-dimensional (2d) 
Ising modeLjn which case detailed studies, both (field- 
)theoreticalaa and numericala cl, have been performed to 
clarify the critical properties of the dilute model. By now, 
according to general view, the dilution is considered as a 
marginally irrelevant perturbation, thus the critical sin- 
gularities in the dilute Ising model are characterized by 
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the power laws of the perfect model modified by logarith- 
mic corrections. There is, however, another view which 
interprets numerical data as .giving evidence for dilution 
dependent critical exponentsQ. 

In this paper we try to decide between the conflicting 
views in an indirect Wjay by studying the local critical be- 
havior at a defect lineEl in the dilute model. A defect line, 
which could be located at grain boundaries in real sys- 
tems, represents a marginal perturbation in the 2d pe: 
feet Ising model. According to Bariev's exact solutio: 
the critical exponent Pd, defined via the temperature de- 
pendence of the defect or local magnetization rrid, 
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t={T,-T)/T, 
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is a continuous function of the strength of the defect cou- 
pling Jd. 
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(a) (b) 
FIG. 1. Ladder (a) and chain (b) type defects. 
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For a chain defect, see Fig. lb, one gets 
2 

/3d = ^ arctan^ Kc, Kc = exp [-2( - J)/Tc] , (2) 

whereas for a ladder defect, see Fig. la, Pd is given by 

. 2 2 tanh(J/r,) 

Pd = ^ arctan ki, m ^ - — , , , , (3) 
TT-^ tanh(J<j/Tc) 

where J is the coupHng in the isotropic Ising model. 
We note that the above formulae could be generalized 
to non-isotropic models, as wellEI. As shown recently by 
Pleimling and Selket3, the edge magnetization of three- 
dimensional Ising magnets at the surface transition has 
a similar non-universal critical behavior, which, indeed, 
can be related to the local critical behavior at a defect 
line in the two-dimensional Ising model. 

The exact results on the local critical behavior of the 
2d Ising model in eqs. (||) and (||) are in complete agree- 
ment with a stability analysis of the fixed point ofpthe 
homogeneous system in the presence of a defect lineM'El. 
Under a small perturbation this fixed point is unstable, 
if the critical exponent of the bulk correlation length, i', 
is < 1. Furthermore, a ladder defect with small local 
couplings behaves like two weakly coupled surfaces, and 
ordinary surface critical behavior will result, provided the 
corresponding surface fixed point remaiss, stable against 
a weak coupling between the surfacesOH. This stabil- 
ity condition can be expressed in terms of the surface 
susceptibility exponent of the homogeneous model as 
7i,i < OeI. Applying hyperscahng, one obtains for d = 2, 
7i^i = v — 2/3i, where f3i is the critical exponent of the 
surface magnetization). As one may easily check the cor- 
responding two marginality conditions, for ladder defects, 

J/ = 1 and 7i,i = , (4) 

are both satisfied for the 2d Ising model; the marginality 
is manifested by the defect coupling dependent critical 
exponent in eq. (^). 

In the following, we are going to utilize the above ob- 
servations and study the local critical behavior at defect 
lines in the dilute Ising model. We consider strongly di- 
luted systems, so that the bulk critical region is clearly 
controlled by the random fixed point, and insert the 
line defects as local perturbation. Then the relevance- 
irrelevance criterion for the local critical behavior is ex- 
pected to have the same form as described above, with 
the exponents, and 71,1, referring now to the dilute 
model. Determining the local magnetization exponent 
(3d, at the defect, one may imagine two scenarios: i) Pd 
showing a continuous variation with the defect coupling 
Jd, or ii) Pd staying constant in, at least, some extended 
range of Jd- In the first case, there would be evidence 
that the marginality conditions, see eq. (^), remain valid 
for the dilute model. Otherwise, one might infer that 
the critical exponents v and 71.1 for the pure and dilute 
models are different. 



In what follows we consider a random-bond nearest 
neighbor Ising model on a square lattice where the ran- 
dom ferromagnetic couplings, Ji and J2|T-9Ccur with equal 
probability. That model is self-dua£3, and the self- 
duality point 

tanh(Ji/r,) = exp(-2J2/Tc) , (5) 

corresponds to the critical point, if there is one phase 
transition in the system. Indeed, this assurnption is 
strongly supported by numerical calculations^. In this 
dilute model, ladder and chain defects are then intro- 
duced, where the defect couplings Jd are uniform and 
ferromagnetic. 

To calculate the local critical properties we did exten- 
sive Monte Carlrt-^MC) simulations using Wolff's clus- 
ter flip algorithmt^l. We have considered square lattices 
with L columns and L rows; Jd couples neighboring spins 
in the center column for chain defects, whereas for the 
ladder defect Jd connects spins between the two cen- 
ter columns. Typically we took L = 256, applying full 
periodic boundary conditions and generating about 10''^ 
clusters per realization. The results are then averaged 
over hundred realizations. The statistical errors during a 
MC run in a given sample turned out to be significantly 
smaller than those arising from the ensemble averaging. 
We mention that similar parameters were used in the 
previous study on. the surface critical behavior of the di- 
lute Ising modelliJ, which corresponds to the case with a 
ladder defect, where Jd vanishes. 

In the MC simulations we calculated the average mag- 
netization per column, m{i) = il'^ ^ijl) / L, where the 
sum runs over j = 1, 2, . . . , L. The defect magnetiza- 
tion is then given by rrid = m{L/2). The simulations 
were performed at three values of the dilution parameter 
r = J1/J2 = 1, 1/4 and 1/10 and for several values of 
the defect coupling in the region < Jd/ J2 < 4. 

The magnetization profile m{i) displays at the de- 
fect either a maximum or a minimum depending on the 
strength of the defect couplings, as illustrated in Fig. 2 
for ladder defects. Far from the defect, there is a plateau 
in the profile with the height signaling the bulk magne- 
tization, rrib. The size of the defect region, Id, where the 
magnetization differs substantially from its bulk value, is 
related to the bulk correlation length of the system. 

In the thermodynamic limit, L 00, as the crit- 
ical temperature Tc, see eq. (^), is approached, the 
magnetization profile 111(1) goes to zero as a power-law 
m{i) - tl^'^'\ where P{L/2) = (3d and = (3 for 

\L/2 — i\ > Id, where /? is the usual bulk critical expo- 
nent. To estimate the values of these critical exponents 
from simulational data, one may define temperature de- 
pendent effective exponents 

/?(z)e//=dln[m(i)]/dln[t] , (6) 

which are approximated by using data at discrete tem- 
peratures, say, t + At/2 and t — At/2. In the limit of 
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sufficiently small At and t, the effective exponents ap- 
proach the true critical exponents, presuming that the 
system is large enough so that finite-size effects play no 
role. 
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FIG. 2. Magnetization profiles in the dilute Ising model 
with a ladder defect at different couplings Jd/Ji =4.0, 0.625 
and 0.01 at the dilution r = Ji/J2 = l/4 and reduced critical 
temperature t — 0.06. MC systems of size 256 x 256 were 
simulated. 



To avoid finite-size effects, L should be much larger than 
the correlation lengths in the bulk and at the defect line. 
Actually we approached the critical point by calculat- 
ing f3{i)eff for t = 0.15, 0.13, 0.11, 0.09 and 0.07, with 
At ~ 0.02, and then did a linear extrapolation to t = 0. 
The error caused by the extrapolation seems to be rather 
small. Further technical details can be found in Ref. 14. 

Before presenting our findings about the defect line 
problem in the dilute model, we will first consider the 
perfect model with random defect couplings. The aim of 
this part of our investigation is to clarify, whether ran- 
dom defects could lead to varying local exponents. The 
two random, ferromagnetic couplings in the defect line, 
J and Jd, are assumed to occur with equal probability, 
where J is also the coupling in the rest of the system. 
Results about the local magnetization exponent Pd for 
various values of the ratio Jd/J are shown in Fig. 3, 
both for chain and ladder defects. The error bars in Fig. 
3 take into account the sample averaging and the extrap- 
olation. Only a few typical error bars are shown. 




FIG. 3. Local magnetization exponent, f3d, of the pure Ising 
model with randomly distributed couplings in the chain and 
ladder type defects, see text. The lines are guides to the eye. 

In Fig. 4, results about the critical exponent of the lo- 
cal magnetization Pd for the dilute model, with random 
couplings Ji and J2, and the defect line, with uniform 
coupling Jrf, are depicted. 




FIG. 4. Local magnetization exponent, Pd, of the random 
Ising model at different dilutions (r = J1/J2 = 1 diamond, 
r = 1/4 square, r = 1/10 circle) for ladder defects. Findings 
on the chain defect are shown in the inset. The exact results 
for the pure system with r = 1 are denoted by full line; the 
broken and dotted lines are guides to the eye. 



For Jd/J — 1, the critical exponent of the perfect 
model, 13= 1/8, is reproduced quite accurately. For other 
values of that ratio, one observes a non-universal critical 
behavior with (3d varying continuously with Jdj J m ac- 
cordance with the above marginality conditions.- It may 
be interesting to note the touching of the curves for the 
ladder and the chain defects at = J . 



In the case r = J1/J2 = 1, our data are in good agree- 
ment with the exact results, see eqs. (^) and (^) for chain 
and ladder defects. In the dilute case, Ji ^ Ji, Pd is seen 
to vary continuously with the strength of the defect cou- 
pling Jd- For a fixed value of Jd, the defect energy-density 
increases, relative to the average bulk value, for decreas- 
ing value of r = Ji/ Ji- Therefore, generally, there is an 
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increasing local order at the defect, which is connected 
to a decreasing value of the defect exponent, (3d- This 
argument, however, seems to be not valid for the ladder 
defect with Jd ^ J2- In this limit one has effectively a 
chain defect with random couplings Ji + J2 , with proba- 
bility 1/2, as well as 2Ji and 2J2, each with probability 
1/4. Then, shown in Fig. 4, (id{Jd) is increasjiiig-.with 
increasing dilution.- For = 0, one recoversoEj the 
surface critical exponent f3d = Pi = 1/2. 

Another limiting situation is obtained for a chain de- 
fect with zero defect bond, = 0. Then the problem 
is equivalent to a ladder defect with three random cou- 
plings, which depend on Ji and J2- As seen in the inset 
of Fig. 4, such random couplings could also lead to a 
non-universal behavior. This observation is in agreement 
with our findings on the perfect model in Fig. 3. 

To summarize, we considered uniform ladder and chain 
defects in two-dimensional dilute Ising models and deter- 
mined the critical exponent of the defect magnetization. 
The exponent was found to be a continuous function of 
the defect coupling. Assuming that the first stability cri- 
terion mentioned above holds for the dilute case as well, 
one gets for the critical exponent v of the bulk random 
Ising model the borderline value v = 1. Accordingly,|jpne 
could rule out > 1, as had been suggested beforeQ in 
the context of dilution dependent bulk critical exponents. 

In conclusion, we suggest that the non-universal crit- 
ical behavior is related to the borderline values of the 
critical exponents of the bulk dilute model, as given in 
eq. (||). Consequently, one obtains v = \ and 71^1 = 
(implying /3i = 1/2, in agreement with Ref. 14), both for 
the perfect and dilute two dimensional Ising models . 
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